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Abstract 

Random Wavelet Series form a class of random processes with multi- 
fractal properties. We give three applications of this construction. First, 
we synthesize a random function having any given spectrum of singulari- 
ties satisfying some conditions (but including non-concave spectra). Sec- 
ond, these processes provide examples where the multifractal spectrum 
coincides with the spectrum of large deviations, and we show how to re- 
cover it numerically. Finally, particular cases of these processes satisfy a 
generalized selfsimilarity relation proposed in the theory of fully developed 
turbulence. 



1 Introduction: Random Wavelet Series 

The story of random Fourier series largely coincides with the development of 
harmonic analysis in the twentieth century, starting with the pioneering work 
of Borel, later developed by Wiener, Salem and Zygmund, to quote but a few, 
and culminating with the famous book of Kahane [?] which made this theory 
accessible to a large audience. Since the mid 80's, wavelet bases have proved to 
be a preponderant alternative option to the trigonometric system, in order to 
analyze and synthesize functions and signals. Therefore, a very natural problem 
is to study random wavelet series. Surprisingly, it turns out that such series have 
properties that differ widely from those of random Fourier series. This difference 
is particularly striking when one considers pointwise Holder regularity. Recall 
that, if a > 0, a function /, defined on K is C a (xo) if there exists a polynomial 
P of degree at most [a] such that 

\f(x)-P(x-x )\<C\x-x \ a . (1) 
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The Holder exponent of f aX xq, denoted by hf(xo), is defined as 

hf(x ) ■= sup{a : / £ C a (x )}- 

Under very general assumptions, random Fourier series with independent coef- 
ficient have everywhere the same Holder exponent. (The differences that appear 
in the modulus of continuity at different points and allow to draw a difference 
between slow points and fast points are logarithmic corrections in the right hand 
side of QJ, of the form log(|x — xol) 7 ; compare for instance in [?] the results 
concerning the uniform modulus of continuity given by Theorem 2 in Chapter 7 
with the results concerning irregularity everywhere in Section 6 of Chapter 8.) 

By contrast, random wavelet series with independent coefficients have a 
Holder exponent which is a highly irregular random function: Let 

Eh := {x, f has Holder exponent h at a;} ; 

if the distributions of the wavelet coefficients depend only on the scale of the 
wavelet (and some mild additional hypotheses) , the sets Eh are non-empty when 
h takes values in an interval of non-empty interior [h m i ni h max ], in which case 
the Eh are random fractal sets, see [?]. The Hausdorff dimension of Eh, which 
is denoted by d{h) is called the spectrum of singularities of the sample path (we 
use the traditional convention dirny(0) := — oo). In this paper we will recall 
the previous results concerning random wavelet series, and continue this study 
by showing in particular that the model supplied by random wavelet series is 
compatible with several turbulence models that have been proposed in the past. 

There are many real-life situations where, although the independence con- 
dition is not necessarily satisfied, good estimates or models are known on the 
distribution of wavelet coefficients. Let us mention a few examples: 

• Several authors (Buccigrossi et al. [?], Huang et al. [?], Mallat [?], Si- 
moncelli [?], Vidakovic [?]) have studied the statistics of the wavelet co- 
efficients of large collections of natural images and observed that these 
statistics are highly non-Gaussian. Exponential power distributions (of 
density Ce- A W) fit very well these statistics. 

• Cascade-type models for the evolution of the probability density function 
of the wavelet coefficients through the scales have been proposed to model 
the velocity in the context of fully developed turbulence (these models 
were initially proposed by Castaing et al. [?] for the increments of the 
velocity, and then fitted to the wavelet setting by Arneodo et al. [?]). 
Random multiplicative models have also been considered in statistics, see 
Vidakovic [?]. 

• Bayesian inference techniques based on a priori models for the distribu- 
tions of wavelet coefficients at each scale have been widely studied (see for 
instance Abramovich et al. [?], Johnstone [?], Miiller et al. [?]) to improve 
the usual wavelet-based denoising algorithms by using the additional in- 
formation supplied by the distributions of wavelet coefficients (instead of 
using only the Besov regularity for instance). 
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• In multifractal analysis, several formulas (Levy-Vehel et al. [?], Evertsz et 
al. [?], Meneveau et al. [?], Riedi [?]) were also proposed in order to derive 
spectra of singularities from distributions of increments of the function. 
These formulas are referred to as large deviation multifractal formalisms; 
they can be easily extended to a wavelet setting. 

We will explore how random wavelet series can play a role in some of these 
models. Let us first recall their definition and main properties. We consider 
functions on T := R/Z (1-periodic functions). Let ip be a mother wavelet such 
that the periodized wavelet family 



form, together with the constant function x i— > 1, an orthogonal basis of L 2 (T) 
(but not orthonormal, the L°° normalization being more convenient for our 
purpose) . The wavelet coefficients of a function / are 



In the following, we will assume that the wavelet ip belongs to the Schwartz 
class, which will simplify the statements of the main results. Note however that 
the regularity results that we will state remain partly valid when using a wavelet 
with limited regularity: In this case, one has to assume that the maximal Holder 
regularity of the process (that will be denoted by h max ) is strictly smaller than 
the uniform Holder regularity of the wavelet. 

1.1 Random wavelet coefficients 

Definition 1. A periodic distribution f is a random wavelet series (RWS) if 
its wavelet coefficients Cj t k in the basis above satisfy the following requirements: 

1. Vj, the Cj,k (k G {0, . . . , 2° ' — 1}) are identically distributed random vari- 
ables; the probability distribution of — log2 ^j,k\) j s denoted by p^; it is 
defined on R U {+oo}; 

2. the Cj^k (j e N, k E {0, . . . , 2 3 ' — 1}) are independent; 

3. there exists 7 > such that 



is strictly negative for a < 7. 

The third requirement is necessary in order to be sure that the series Y^, Cj,ki ) j,k 
is convergent in the sense of distributions. The function p thus defined is called 





a) := inf limsup 

£>0 j^ + oc 



log 2 (2 3 Pj([a-e,Q + e])) 
3 



(2) 
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the upper logarithmic density of the process. It is upper semi-continuous, but 
not necessarily monotonous. We do not make any other assumption on the 
probability measures py, note that p-({+oo}) is simply the probability that 
Cj.k = 0. 

1.2 Histograms of wavelet coefficients 

Consider now an arbitrary function (or distribution) /, which, for instance, can 
be a sample path of a stochastic process. Let Nj(a) := # {k, \Cj,k\ > 2~ aj }, 
and 

i . ... log 2 (N j (a + e)-N j (a-e)) , n , 
p(a) := mf hmsup bzy Jy - J — — . (3) 

Without any additional assumption, we expect p to be a random function, 
whereas, if / is a RWS, p is deterministic. A first result of [?] links these two 
functions. 

Theorem 1. Let f be a RWS, and 

W := < a,Ve > 0, ^ 2 j p j ([a - e, a + e]) = +oo > . 

[ j£N J 

With probability one, for all a, 



p(a) 



p(a) if a G W: 
— oo eZse. 



One immediately checks that p(a) > ^> a E W and p(a) < => a £ W. 
If p(a) = 0, then p{a) can be either (almost surely) or — oo. 

1.3 Spectrum of singularities 

A second result gives the spectrum of singularities for a RWS Naturally \{W = 0, 
then / is almost surely globally C°° . Otherwise, let us define 

/imin := inf(W), 
and, assuming that 3a, p(a) > 0, 

, ( P(a) 

ftmax := SUp 



*>0 a 

Note that, because p is upper semi-continuous, p(h min ) > 0. 
Theorem 2. Let f be a RWS Almost- surely, 

• the almost- everywhere Holder exponent is h max ; 
for all h £ [/i min , hmax], d(h) = hsup ae(o h] 
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• d(h) = — oo else. 

See Figure for an illustration. 

In [?] it was implicitly assumed that there exists a such that p(a) > 0. 
However, it may happen that this condition does not hold, and yet W is not 
empty. Consider for instance the RWS defined by fixing ag > 0, and Vj G 
N, Pj({a }) = j2~i and Pj({+oo}) = 1 - j2~i (here W = {a }). In that 
(degenerate) case, we get an almost sure "flat" spectrum. 

Proposition 1.1. Let f be a RWS such that Va, p{a) < 0, but W ^ 0. Define 
hrain as above. Then, almost surely, 

• the almost- everywhere Holder exponent is +00; 

• for all h > h mln , d(h) = 0; 

• d(h) = —00 else. 

Proof. First, according to Theorem^ p(l) < implies a. s. ,0(7) = — 00, which 
means that with at most finitely many exceptions, \Cj^\ < 2~ 7J , 
For a > let K j (a) := {k, \C j)k \ > 2"^'} and if d'< 1, let 

E j (a, d) := [j {k2~ j - 2"*', k2~ 3 + 2~ dl ), 

and 

E(a,d) := limswp E J (a, d). 

Note that E(a,d) is increasing in a and decreasing in d. Because p(a) < 
0, Vd > 0, by the Borel-Cantelli lemma, almost surely E(a, d) has Lebesgue 
measure C{E{a,d)) = 0. Then, with 

m>l 

almost surely C(E) =0. If x E, then for all m > 1, for all j, k with at most 
a finite number of exceptions, either \Cj.k\ < 2~ m i or | a? — fc2 | > 2~^, in 
which case \Cj,k\ <\x — k2~ j \ m . Using the classical wavelet characterization of 
pointwise regularity, this proves that / <E C°°(x). 

The last two points are similar to Theorem (in the case where p(a) > 
happens only for some a > ho > h mm ). 

□ 

To conclude this overview, let us mention the following result, also proved in 
[?] . A function / is called a uniform Holder function if there exists e > such 
that / e C* e (T). 

Proposition 1.2. The spectrum of singularities of any uniform Holder function 
f satisfies the inequality 

d{h) < h sup ffl. (4) 

ctS(0,/i] a 
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2 Synthesis of multifractal processes 



When a multifractal model is proposed, a natural step is to construct a function 
or a random process having a given spectrum of singularities d(h). According 
to Theorem [3 the candidates for RWS spectra are right-continuous functions 
satisfying the following properties: d < 1; d > on an interval [h m in,h m ax] 
where the function h t— > ^f^- is increasing; outside this interval d(h) — — oo; and 
d{h max ) — 1. These conditions are also sufficient. 

Proposition 2.1. Let d{h) satisfy the conditions listed above. Take, for all 
j G N. for all a G K, 

Pj{da) := iM?) 2 i(^)-i) dQ 

"'max 

and p^d+oo}) := 1 — J ' lmax Pj(da). Then Pj is a probability measure on RU 
{+oo}, and its upper logarithmic density is p{a) = d(a). 

Proof. Note that d(a) < It follows that 

hn 



J^p j (da)=J Pj(da) 



<^2-i[ h ™2^da 



< 1, 



which ensures, together with the definition of Pj({+oo}), that p - is a probability 
measure on M. U {+00} . 

Let us now compute the upper logarithmic density. We have 

^ '- + d a < " : < — + d(a + e), 

3 3 3 

\o g2 h £^ P] {t)dt) 
d(a) < limsup : — < d(a + e) 

i^+oc 3 

hence, letting e — > and using right-continuity, p(a) = d(a). □ 

The almost sure spectrum of singularities of the corresponding RWSis then 
d. To synthesize it, for all j G N we draw independently V random variables ctjk 
with law pj (using the rejection method if necessary) , and let Cj t k := Xjk^~ Joljk 
(Xjk is an arbitrary sign or phase). Then 

= EE C jtk ^ h {x) (5) 
jeN fc=o 

is a realization of the process. 
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This allows in particular to show processes with non-concave spectra. An 
example is given on Figure^with d(h) := (h— ^) 2 on [5, |] and — 00 elsewhere. 
The wavelet used is Daubechies 10 (extremal phase). 

0.008 I 1 1 1 1 1 1 1 



0.006 



0.004 



0.002 







-0.002 



-0.004 




-0.006 1 ' 1 ' ' 1 ' 1 

5000 10000 15000 20000 25000 30000 35000 

Figure 1: A multifractal function with non-concave spectrum 



3 Spectrum of large deviation 

For a function given in a sampled form (signals, images), one often wishes to 
compute its spectrum of singularity. Just applying the definition (identifying 
the sets of common Holder exponent and computing their Hausdorff dimension) 
is clearly not feasible, so numerical procedures called multifractal formalisms 
have been proposed instead. In general, a multifractal formalism does not yield 
the correct spectrum, but an upper bound, and is proved valid only on certain 
functions. For instance, the "classical" multifractal formalism, derived from 
the original ideas of Parisi and Frisch [?], consists in computing the so-called 
structure function 

j^ + oo —j 

and taking its Legendre transform 

di(h):=in£hq-T(q) (6) 

9 

as an estimation for d{h). It was proved by Jaffard [?] that for any function 
/, if q c is the (only) solution to r(g c ) = 0, and if the infimum in © is taken 
for q> q c , then d(h) < di(h); moreover, equality in (jjjj holds for selfsimilar 
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functions (Jaffard [?]). Note, however, that this cannot be true for all func- 
tions; in particular the right-hand term of is concave, whereas in general the 
spectrum is not. 

One can also prove that, for any function (even for a tempered distribution), 
for all h, 

h sup ^ < inf hq-T(q); (7) 
this implies that 10} is sharper than 0. Actually, 

d 2 {h) :=h sup ^ 

ae(0,h] a 

is equivalent to the spectrum of large deviations of [?] ; we know from Theorem[2 
that for RWS it coincides with the spectrum of singularities. 

Because of 0, this new multifractal formalism will be valid whenever the 
classical one is valid; moreover, since the spectrum thus obtained is not necessary 
concave, its domain of validity is strictly larger than the classical one. Note that 
p{a) may not be easy to compute numerically, because © involves a double 
limit. But if we define 

A(a) := hmsup . ' , 

which is increasing, we can show that its upper closure A (whose hypograph 
is the closure of the hypograph of A) satisfies A (a) = sup Q / <Q p{a'). Hence 

d2(h) — sup aG( -Q h \ as well, which is easier to compute. 

We tested this algorithm on the process that we synthesized in § with 
the same parameters as on Figure except that it was computed with 2 22 
points to get a sufficient scale range. The analyzing wavelet is Daubechies 3 
(extremal phase), different from the synthesizing wavelet. Implementation is 
straightforward; we used a simple linear regression on the 10 largest scales to 
compute X(a), and then da{h). Results are shown on Figure^] 



4 Generalized selfsimilarity 

The model for fully developed turbulence proposed in Castaing et al. [?] asserts 
that the velocity field is a random process X, with increments at scale I following 
a law of density Pi, and that if I < L, 

Pi(x) = J g lL (u)e- u P L {e- u x)du, (8) 

where the selfsimilarity kernel satisfies giL — giv * gi'L for I < V < L. This is a 
generalization of the notion of selfsimilar process, because taking giz := S Hln i 
in ® yields 



Pi(x) = 
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Figure 2: Computed spectrum of singularities: numerical (+) and theoretical 
( — ) results. 



The construction of such a process for general g is still an open prob- 
lem, but a discrete approach can be done using wavelets. Assuming that for 
all k, Cj.k=X 2 ~j , and that, at j fixed, the common probability density for 
— log 2 (|Cj,fc|) is pj, (SJ becomes for j > J 

Pj = g jJ * Pj> (9) 

where Gjj{u) :— .92-32- ' ( — u ) must satisfy Gjj — Gjj' *Gj'j for j > j' > J. One 
can furthermore assume that Gjj depends only on j — J, in which case Gjj — 
Q*(j-J)_ Remark that, with the notation of Definition Pj(da) — jp~j(ja)da. 
A RWS satisfying © can be obtained by simply taking Pj(da) = jGjo(ja)da. 
Then its spectrum of singularities can be computed almost surely and, thanks 
to Theorem it satisfies the multifractal formalism given by © , where the 
infimum is taken for q > q c . 



Example 1: G = A r m cr 2, with to > a log ^ . For all a, 

(a — to) 2 



p(a) = 1 -log 2 (e) 



2(7 2 



The wavelet coefficients follow a log-normal law; the spectrum of singularities 
is a segment of a parabola followed by a segment of a line. 
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Example 2: G = S aa * 7^/3, with v, (3 > 0; ao > a*{v, (3), which is the largest 
solution to 1 + ^ log 2 (— a*) + /31og 2 (e)a* + ^log 2 (^f ) =0. For a > a , 



p(a) = 1 + ^log 2 (a - a ) - /31og 2 (e)(a - a ) + ^log 2 




Example 3: G = S aa * p c , where p c is a Poisson distribution with parameter 

c and ao > a*(c), which is the solution to 1 — clog 2 (e) — a* log 2 (t^t) = 0. 

This kernel was proposed by Dubrulle [?] and She and Waymire [?] in the study 
of fully developed turbulence. For a > ao, 




Figure 3: Upper logarithmic density (thin curve) and spectrum of singularities 
(bold curve) for the Poisson kernel (c = 1, ao = 0). 



It is remarkable that such selfsimilarity kernels can be obtained in the frame- 
work of RWS, since they are usually expected to be the signature of cascade 
models on the wavelet coefficients which display strong correlations between 
these coefficients (see [?]). Our examples show that it is not the case: Correla- 
tions between wavelet coefficients cannot be inferred from the particular shape 
of the p.d.f. of the wavelet coefficients at each scale. Note however that one 
possible option, in order to derive some information on these correlations, is to 
study how statistics of local suprema of the wavelet coefficients behave, see [?]. 
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